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Abstract 

As a powerful tool in scientific computation, Mathematica offers us algebraic com- 
putation, but it does not provide functions to directly calculate commutators in 
quantum mechanics. Different from present software packets to deal with noncom- 
mutative algebra, such as NCAlgebra and NCComAlgebra, one simple method of 
calculating the commutator in quantum mechanics is put forward and is demon- 
strated by an example calculating SO (4) dynamical symmetry in 3 dimensions 
Coulomb potential. This method does not need to develop software packets but 
rather to directly write program in Mathematica. It is based on the connection be- 
tween commutator in quantum mechanics and Poisson bracket in classical mechanics 
to perform calculations. Both the length and the running time of this example are 
very short, which demonstrates that this method is simple and effective in scientific 
research. Moreover, this method is used to calculate any commutator in Hamilton 
system in principle. In the end some deficiencies and applications are discussed. 
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1 INTRODUCTION 



Today scientific computations are classified into two kinds. One is pure nu- 
merical computation, the other is algebraic computation. The latter treats 
with symbolic computation rather than simple numerical computation, thus 
it holds more general applications in scientific research. Among much software 
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to implement symbolic computation, Mathematica, developed by Wolfram Re- 
search Inc., is one of the most popular, influential and active software in aca- 
demic and industrial world [1]. It offers us powerful scientific computational 
facilities, algebraic computation, high-quality graphical visualization and user 
interfaces. Moreover, its basic system is developed in C program language, 
and it is easy to be transplanted into other programs and languages. 

Mathematica has been extensively applied in scientific research, especially in 
mathematics, physics and chemistry. For instance, it has been used to solve 
Schoringer equation for bound states [2], calculate tree-level multi-particle 
electroweak production process [3], and perform numeric computation of dyadic 
Green's function [4]. But there still exist many other problems on symbolic 
computations to be simplified. Calculating commutators in quantum mechan- 
ics is one of such problems. 

The emergence of the commutator is a fantastic event in physics [5]. In classical 
mechanics there does not exist noncommutative quantity, while in quantum 
mechanics, there are two kinds of quantities: Q quantity and C quantity. Q 
quantity denotes noncommutative quantity. For example, when two quantities 
A, B do not meet commutative law, i.e. [A, B] = AB — BA ^ 0, we may say 
that A, B are two Q quantities. Otherwise, we may say that A, B are two 
C quantities. In quantum theory, especially in quantum mechanics, quantum 
group and statistical physics, we often have to calculate many commutators. 
Today the research of the commutator has initiated noncommutative algebra 
[6] and noncommutative geometry [7,8], which are all hot topics in physics. 

However, in practice, to perform these symbolic calculations is a tedious and 
mechanical job for physicists, so it is necessary to seek a simple and effective 
method of calculating operators with computer. It should be noted that even 
though Mathematica does provide symbolic computations, it does not include 
calculating commutators yet. Therefore, a direct idea is to develop source 
codes in Mathematica to compute commutator or anti-commutator. Along 
this routine, scientists have developed some software packets to deal with 
noncommutative algebra [9-11], such as NCAlgebra [9] and NCComAlgebra 
[10], to perform symbolic computations. On one hand, to develop, install and 
use these professional packets need many skills and knowledge of computer, 
which might be hard work for most physicists. On the other hand, in these 
packets, we should set the basic commutative relations firstly, and then by us- 
ing computer, we simplify these operators. The process makes it inconvenient 
to calculate commutator with this software. Moreover, some software is com- 
mercial and it would not be used free. Hence, if one simple way to calculate 
the commutator is put forward, it will bring us convenience in practical re- 
search. After all, we need a simple method. Actually, noted that Mathematica 
provides differential calculation, and with the relation between commutator 
in quantum mechanics and Poisson bracket in classical mechanics, it might be 
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used to calculate the commutators. Based on this idea, one simple, easy new 
method of calculating commutators in Hamilton system is presented and its 
process is illustrated by one example in this article. The method does not need 
to develop any software packet or modify any source code in Mathematica and 
it executes calculations only with functions in Mathematica software itself. It 
should be noted that this method is very simple and it does not need other 
extra skills or knowledge of computer. Moreover, this method can be extended 
into solving more problems. In the end some defects and applications of this 
method are discussed. This method is expected to simplify many practical 
computations in further scientific research. 



2 POISSON BRACKETS AND COMMUTATORS 



Now let us recall the relation between Poisson bracket and commutator 

{A,B} = ±[AB], (1) 



where [ ] donates commutator, {} donates Poisson bracket [12], i is the imagi- 
nary number unit, and h is Planck's constant. This relation gives us the con- 
nection between the quantum mechanics and classical mechanics. Assume a 
Hamilton system H(pi, where pi, qi are canonical coordinates and their con- 
jugate momenta respectively. For two mechanical quantities A(pi, B(p i7 
the Poisson bracket is defined by [13] 

^ f 8AdB OA OB 
{AB} = U^ wrWiWi ). (2) 



From Eq.(l) and (2), we will obtain 



[A,B]=ih'£ 



dAdB dAdB 

dq { dpi dpi dqt 



(3) 



Eq. (3) provides some clues to calculate commutators. It should be noted that 
the left side of eq.(3) is noncommutative algebra, but the right-hand side 
becomes a polynomial in classical mechanics and there does not exist non- 
commutative algebra, i.e. eq.(3) give us a tool to change noncommutative al- 
gebra calculations into commutative algebra calculations. When we compute 
commutators in quantum mechanics, we firstly rewrite these commutators 
into Poisson brackets in classical mechanics. Then we could use the functions 
in Mathematica directly to calculate those Poisson brackets, and at last we 
rewrite the results into quantum operators and get calculating results. Because 
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Mathematica provides powerful functions, such as Factor, Simplify, D[f,x] [1], 
we can effectively finish complex calculations in a few minutes. In the next sec- 
tion we will describe this method by computing SO (4) group in 3 dimensions 
Coulomb potential. 



3 CALCULATING SO(4)GROUP IN 3 DIMENSIONS COULOMB 
POTENTIAL WITH MATHEMATICA 



Assume a particle in Coulomb potential V(r) = — -, we try to calculate the 
dynamical symmetry group SO (4). In 3 dimensions Coulomb potential, there 
exist three conserved quantities, Hamiltonian H, angular momentum L and 
Runge-Lenz vector R — |(P x L — L x P) — e r [14]. In order to know whether 
angular momentum L and Runge-Lenz R construct SO (4) Lie group, we should 
calculate the commutators [H, R], [R, R], [R, L], [L, L\. Firstly we calculate the 
commutator [R,R\. In fact from the symmetry, we might only calculate the 
commutator [R x ,R y ], where R x ,R y are Runge-Lenz vector R along x and y 
axis respectively. Generally, calculating these commutators is a tedious job, 
but we will show that with Mathematica it becomes a piece of cake. 

First. Simplify the commutator to be calculated with the algebra relation 
PxL + LxP = 2ihP. Then, we obtain 

[R x , R y ] — [{P x L — ihP) x — e x ,(P x L — ihP) y - e y \. (4) 



In most circumstances, we should simplify the commutators before we calcu- 
late these commutators. 

Second. Using eq.(l) and eq.(4), we rewrite above commutator into Poisson 
bracket as follows: 

[Rx, R y ] = ih{(P x L — ihP) x — e x ,(P x L — ih P) y - e y }, (5) 



here we should note the definitions (PxL) x = ih(l z p y — l y p z ), l x = ih(yp z — zp y ) 
and e x = r , e v = r etc. 

Third, eq.(5) also can be rewritten into 



o 

[R x , R y ] = ih E(^(( p x L - ihp )* - e ^dp-tt p xL ~ lhp )y ~ e v) ( 6 ) 



-^-((P XL- iHP) x - e x )ji(P xL- ihP) y - e y )), 
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where i = x,y,z; pi = p x ,p y .p z . Using the function in Mathematica(see Ap- 
pendix), we could directly calculate the right-hand of eq.(6) as follows 

2 

[R x , R y ] = -ih(p y x - p x y)(P 2 - -). (7) 

Similarly, we could calculate other commutators with Mathematica. 

[H,R x ] = 0, (8) 

z 

[L x , R y ] = -ih(p x p z x + p y p z y - zp 2 x - zp\ - ihp z - -), (9) 

[L x , L y ] = -ih(xp y - yp x ). (10) 

Fourth, rewrite above results into the standard quantum mechanics operators 
and simplify these results. We get 

[R x ,Ry] = -2ihHL z , (11) 

[H,R] = 0, (12) 

[L x ,R y ]=ihR z , (13) 

[L x ,L y ]=ihL z . (14) 

Up to now we have finished calculating such tedious commutators. From above 
relations, we can see that the Rungle-Lenz vector R and angular momentum 
L construct the SO (4) dynamical symmetry, which is higher than geometry 
symmetry SO(3). This is a well-known result in quantum mechanics [15,16]. 
What we want to do here is not to prove or check this result, but to demon- 
strate that our method is efficient. It should be noted that the program is very 
short(only 20 lines) and the running time is only half a minute (PC 300, 64MB 
RAM)! In contrast to calculating above commutators by hand, it is faster and 
more exact. The program and some expeditions are listed in Appendix and 
there the function Commutator3D can be used to compute any commutator 
in 3 dimensions Hamilton system. 



4 DISCUSSIONS 



Because the connection between the commutator and Poisson bracket only 
helps to calculate the Hamilton system in quantum mechanics, our method 
might be unsuitable for the second quantize system. However, it does pro- 
vide a simple and powerful method of calculating or checking commutators 
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on dynamics symmetry and Lie algebra in quantum mechanics and classical 
mechanics, especially when the Hamiltonian contains complicated potentials. 
In principle, this method is used to calculate any commutator in Hamilton 
system. We test our method in several problems, such as calculating Lie al- 
gebra or Lie group in isotropic harmonic potential [17], screened Coulomb 
potential and screened isotropic harmonic potential [18]. All of the results 
have demonstrated that this method is simple, effective and robust. We ex- 
pect this method could help physicists solve such problems easily. Moreover, 
this method is easy to be applied into other mathematical software, such as 
Maple, Matlab. Meanwhile our method indeed shows that for calculating most 
commutators in Hamilton system, Mathematica could solve these problems by 
itself. 

Of course, it is interesting and valuable to develop these software packets [9- 
11] to perform symbolic calculations, because they will release us from the 
bondage of mechanical and tedious work, and help us conduct real scientific 
research. In contrast to using these software, our method seems to be simpler 
and easier for most physicists. Of course, this method could not take the 
place of these software packets and itself is only supplement to these software 
packets. Actually we also hope that our method could give some clues in 
developing noncommutative algebra software and make it more powerful. 
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APPENDIX 



(*This program is used to calculate the SO(4) group in 3D Coulomb potential 
V(r) = -1/r *) 

(* Comumutator3D is a function to calculate commutators in 3 dimensions, 
where f and g are variables to be computed *) 

(* pi and qi (i=l,2,3) are canonical coordinates and their conjugate momenta 
respectively *) 

(*ham is Hamiltonian of system, h is Planck's constant, 1 is angular momen- 
tum and r is Rungz-Lenz vector*) 

Commutator3D[f_, g_]:= i h {D[f, qi ] D[g, Pl ]- D[f, Pl ] D[g, qi ]+D[f, q 2 ] D[g, 
pa] -D[f, p 2 ] D[g,q 2 ]+ D[f, q 3 ] D[g, p 3 ]- D[f, p 3 ] D[g, q 3 ] } 

ham={(p 1 "2+p 2 "2+p 3 "2)/2}-l/Sqrt[qr2+q 2 "2+q 3 "2)] 
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ii=q2P3-q3P2 
l2=q3Pi-qiP3 
l3=qiP2-q2Pi 



ri=l3P2-l2P3-i h pi-qi/Sqrt[qi"2+q 2 "2+q 3 "2)] 
r2=liP3-l3Pi-i h p 2 -q 2 /Sqrt[qi"2+q 2 " 2+q 3 "2)] 
r 3 =l 2 pi-lip 2 -i h p 3 -q 3 /Sqrt[qi"2+q 2 "2+q 3 "2)j 

Si =Commutator3D [r i ,r 2 ] 

s 2 = Commut at or 3D [r i , ham] 

s 3 =Commutator3D[li,r 2 ] 

S4=Commutator3D[l 1 ,l 2 ] 

fi=Simply[si] 

f 2 =Simply[s 2 ] 

f 3 =Simply[s 3 ] 

f4=Simply[s4] 
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